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ABSTRACT 

We present the resnlts of three-dimensional global resistive magnetohydro¬ 
dynamic (MHD) simnlations of black hole accretion flows. General relativistic 
effects are simnlated by nsing the psendo-Newtonian potential. Initial state is 
an eqnilibrinm model of a torns threaded by weak toroidal magnetic helds. As 
the magnetorotational instability (MRI) grows in the torns, mass accretes to the 
black hole by losing the angnlar momentnm. We fonnd that in the innermost 
plnnging region, non-axisymmetric accretion flow creates bisymmetric spiral mag¬ 
netic helds and cnrrent sheets. Mass accretion along the spiral channel creates one 
armed spiral density distribntion. Since the accreting matter carries in magnetic 
helds which snbseqnently are stretched and amplihed dne to diherential rotation, 
cnrrent density increases inside the channel. Magnetic reconnection taking place 
in the cnrrent sheet prodnces slow mode shock waves which propagate away from 
the reconnection site. Magnetic energy release in the innermost plnnging region 
can be the origin of X-ray shots observed in black hole candidates. Nnmerical 
simnlations reprodnced soft X-ray excess preceding the peak of the shots. X-ray 
hardening at the peak of the shot, and hard X-ray time lags. 

Subject headings: accretion,accretion disks - black hole physics - instability - 
MHD 


1. Introduction 


Magnetic fields play essential roles in various activities associated with accretion disks. 
The most important role is the angular momentum transport which enables the accretion of 
the disk material. Since Balbus & Hawley (1991) pointed out the importance of the magneto- 
rotational instability (MRI) in accretion disks, the nonlinear growth of the instability has 
been studied by local three-dimensional magnetohydrodynamic (MHD) simulations (e.g., 
Hawley, Gammie, & Balbus 1995; Matsumoto & Tajima 1995; Brandenburg et al. 1995) 
and global MHD simulations (e.g., Armitage 1998; Matsumoto 1999; Hawley 2000; Machida, 
Hayashi, & Matsumoto 2000). 

When general relativistic effects are taken into account, accretion flow in the innermost 
region of black hole accretion disks changes from a circularly rotating flow to a radially in- 
falling flow. There exist some radius below which almost all energy is advected inward into 
the black hole. The radial flow is transonic around the last stable orbit. Steady models of 
black hole accretion disks including this transonic region were constructed under the pre¬ 
scription of a-viscosity for optically thick disks (Muchotrzeb & Paczyhski 1983; Matsumoto 
et al. 1984; Abramowicz et al. 1988) and for optically thin disks (e.g., Matsumoto, Kato 
& Fukue 1985; Narayan, Kato & Honma 1997). Global three-dimensional MHD simulations 
of black hole accretion flows have been carried out by Hawley (2000), Hawley (2001), Haw¬ 
ley & Krolik (2001), Haley, Balbus & Stone (2001), Armitage, Reynolds & Ghiang (2001), 
Reynolds & Armitage (2001), Hawley & Krolik (2002), Hawley & Balbus (2002), and Krolik 
& Hawley (2002) by using the pseudo-Newtonian potential (Paczyhski & Witta 1980). In 
black hole accretion flows, they reproduced radial flow prohles similar to those obtained by 
assuming the phenomenological a-viscosity. Hawley & Krolik (2001) showed that the ratio of 
stress to pressure which corresponds to the a parameter in conventional disk models, exhibits 
both systematic gradients and large fluctuations; it rises from 10“^ in the disk midplane at 
large radius to ~ 10 near the midplane well inside the marginally stable radius. Thus the 
efficiency of angular momentum extraction in the plunging region can be larger than that in 
transonic disk models (see also Reylonds & Armitage 2001). 

Another important mechanism in the innermost region of black hole accretion flow is 
the magnetic reconnection which can release magnetic energy even inside the marginally 
stable radius. In differentially rotating disks, magnetic helds are stretched and amplihed 
due to differential rotation. When uniform horizontal held threads an accretion disk, for 
example, magnetic helds twisted by diherential rotation create current sheets which subject 
to magnetic reconnection. Tajima & Gilden (1987) carried out MHD simulations of the 
formation of current sheets and magnetic reconnection in accretion disks. They applied the 
mechanism to quasi-periodic oscillations of dwarf nova. Sano & Inutsuka (2001) carried 
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out three-dimensional local MHD simulations of MRI and showed quasi-periodic release of 
magnetic energy in accretion disks. 

In our previous papers (Machida et ah 2000; Kawaguchi et ah 2000; Machida et 
ah 2001), we showed the developments of turbulent magnetic helds and current sheets 
in accretion disks by three-dimensional global MHD simulations assuming ideal MHD. By 
evaluating the released magnetic energy from the current distribution obtained by ideal MHD 
simulations, Kawaguchi et al. (2000) successfully reproduced 1//" type time variations 
observed in black hole candidates. 

In addition to the fractal-like 1//" time variations, Negoro et al. (1995) pointed out that 
the X-ray time variations from Cyg X-1 show large amplitude X-ray shots. The X-ray shots 
have typical time interval around several seconds and show symmetric time prohle (Negoro 
et al. 1995; Negoro, Kitamoto, & Mineshige 2001; Negoro & Mineshige 2002). Manmoto et 
al. (1996) showed by hydrodynamical simulations based on a-prescription of viscosity that 
inward propagating waves in advection dominated black hole accretion flows are reflected 
around the marginally stable radius and reproduced time symmetric X-ray prohles. They 
attributed the origin of rapid (~ ms) spectral hardening at the peak of X-ray shot to the 
temperature rise due to acoustic wave reflection. 

In this paper, we explicitly include the electric resistivity in basic equations and present 
the results of three-dimensional global resistive MHD simulations which show the magnetic 
energy release following the build up of current sheets in the innermost region of black hole 
accretion flows. 


2. Basic Equations and Simulation Model 

We solved the following resistive MHD equations in cylindrical coordinate system (cu, z) 
by using a modihed Lax-Wendroff scheme with artihcial viscosity. We use the normalization 
c = Tg = 1 where c is the light speed and Vg is the Schwarzschild radius. 


| + VM = o 
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dv 

— + V - VV 
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where p, P, v, B, 0, j, T and S are the density, pressnre, velocity, magnetic helds, gravita¬ 
tional potential, cnrrent density, temperature and specihc entropy, respectively. The specihc 
entropy is expressed as S' = C„ln(P/p''') where is the specihc heat at constant volume 
and 7 is the specihc heat ratio. Joule heating term is included in the energy equation. We 
neglect radiative cooling. We assume the anomalous resistivity rj adopted in solar hare simu¬ 
lations (e.g., Yokoyama & Shibata 1994) which mimics the enhancement of resistivity where 
current density j is large, as follows; 

rj = r]o[max{vd/vc - 1, 0)]^ (5) 


where Vd = j/p, and Vc is the threshold above which anomalous resistivity sets in. We adopt 
Vc = 8 , and rjo = 5 x 10“^. 


The initial condition is an equilibrium model of an axisymmetric MHD torus threaded 
by toroidal magnetic helds surrounding a black hole. We assume pseudo-Newtonian potential 
0 = —GM/ (r — Tg) (Paczyhski & Witta 1980) where G is the gravitational constant and M 
is the mass of the black hole. We neglect the self-gravity of the gas. At the initial state, the 
torus is assumed to have a constant specihc angular momentum L and assumed to obey a 
polytropic equation of state P = Kp'^ where P is a constant. We assume that the torus is 
embedded in hot, non-rotating, isothermal, spherical halo. 


According to Okada, Fukue & Matsumoto (1989), we assume the square of the Alfven 
speed v\ = P^/( 47 rp) = H{pzu‘^y~^ , where is the toroidal magnetic held and P is a 
constant. Using this assumption, we can integrate the equation of motion into a potential 
form; 


z) 


1 

2 (r - 1) 


, 1 2 , 7 2 


Tfe = constant, 


( 6 ) 


where = 'jP/p is the square of the sound speed, r = -(- and Tf, = 4/(^6, 0 ). 

Here, the reference radius zub is dehned as the radius where the rotation speed Ljwb equals 
the Keplerian velocity v-^h = [^b{d(i)/. By using equation ( 6 ), we obtain the 
density distribution as 


max[Tfc + l/(2(r - 1)) - LV(2 zi7 ^), 0] | 

^[ 7 /( 7 - 1)][1 i 


where jSb = is the ratio of gas pressure to magnetic pressure at (zu,z) = 

{zub,0). We set pb = 1. The parameters describing the structure of the MHD torus are 7 , 
Pb, L, and K. In this paper we report the results of simulations for parameters zub = SOr^, 
Pb = 100, 7 = 5/3, L = (cc 7 ^/ 2 )^A . zub/i^b — 1), and K = 0.0005. The density and sound 
speed of the halo at r = zub is taken to be phaio/pt = 10“^ and n^^haio = 3/5, respectively. 
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The number of mesh points are N^, Nz) = (250,64,192). The grid size is Aw = 
Az = 0.1 for 0 < w, 2 : < 10, and otherwise increases with w and 2 :. The outer boundaries at 
w = 150 and at z = 70 are free boundaries where waves can transmit. A periodic boundary 
condition is imposed for the (p-direction. Since we include full azimuthal angle (0 < v? < 27r) 
into the simulation region, we can follow the evolution of low azimuthal wave number non- 
axisymmetric structures. We imposed symmetric boundary condition at the equatorial plane. 
Absorbing boundary condition is applied at the inner boundary by introducing a damping 
rate 

a, = 0.1(1.0-tanh "'~^;^^^"^ ). (8) 

2Aw 

Inside r = 2, the deviation of physical quantity q from initial value go is artihcially reduced 
with damping rate a* as 

= q-ai{q-qo). (9) 

This damping layer serves as the non-reflecting boundary which absorbs accreting mass and 
waves propagating inside r = 2. A small-amplitude, random perturbations are imposed at 
t = 0 for the azimuthal velocity. 


3. Numerical Results 

3.1. Global Structure of Non-Radiative Accretion Disk 

Figure la and hgure lb show the isosurface of the density (p = 0.4) at t = 0 (initial 
condition), and t = 30590 ~ lOto where to = 27iWb/vKb — 3079 is the one orbital time at 
Wb- Owing to the efficient angular momentum redistribution, the torus is deformed into 
a disk and matter accretes to the central part by losing the angular momentum. In the 
outer part of the torus, matter gets angular momentum and expands radially. The angular 
momentum transport is initially driven by the magneto-rotational instability. After the 
non-axisymmetric MRI grows, the disk region becomes turbulent. 

Figure 2a shows the azimuthally averaged density contours and the momentum vectors 
at t = 30590. Circulating motions similar to those observed in ideal MHD simulations (e.g., 
Hawley & Krolik 2001; Machida et ah 2001) appear near the equatorial plane {w 12). 
In the innermost region {w < 5) radial infall becomes signihcant. Near the surface layer 
around ru ~ 7, shock wave front is formed and higher angular momentum gas blows out into 
the disk corona. Figure 2b shows the azimuthally averaged density distribution at t = 0, 
t = 21502 and t = 30590. 

Figure 3a and 3c show the equatorial density distribution. Figure 3b and 3d show 
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magnetic field lines projected onto the equatorial plane. In the outer region (hgure 3a, 3b), 
magnetic held lines are tightly wound and show turbulent structure. In the inner region 
(hgure 3c, 3d), magnetic held lines are less turbulent and globally show bisymmetric spiral 
(BSS) shape. One armed spiral mode dominates in the density distribution. Figure 4 shows 
three-dimensional structure of magnetic held lines at t = 30590 around the surface of the 
innermost region {w < 10) of the disk. White curves show the magnetic held lines. Yellow 
plane shows the equatorial plane and blue surface is the isosurface of the density (p = 0.4). 
Magnetic held lines have signihcant 2 :-components and show helical structure near the surface 
of the disk. Matter swirls into the black hole along these held lines. Bisymmetric spiral 
magnetic helds are created owing to the infall of the disk material. 


3.2. Time Evolution and Angnlar Momentnm Redistribution 

The amplihcation of magnetic helds and its saturation occurs similarly to those reported 
by Hawley & Krolik (2001, 2002). Figure 5a shows the time evolution of the ratio of the gas 
pressure to magnetic pressure (3 = (P) /where the bracket means the volume average. 
The solid curve and the dashed curve show the average in the inner region (4 < cu < 10, 
0 < < 1) and in the outer region (20 < zu < 40, 0 < < 3), respectively. The disk 

stays in a quasi-steady state with /? ~ 10 for time scale much longer than the dynamical 
time scale. Figure 5b shows the time evolution of the ratio of the Maxwell stress to pressure 
as = — (ProP(^/47r)/(P) averaged in the inner region and in the outer region. We found 
aB — 0.1 in the inner region and as — 0.02 in the outer region. Figure 5c shows the time 
variation of the accretion rate at cu = 2.5. Accretion rate increases with time even at the end 
of simulation because much mass is still stored in the original torus {zub ~ 50). It is beyond 
the scope of this paper to continue simulation for time scale long enough for the black hole 
to accrete most of the mass of the initial torus. 

The Maxwell stress in hgure 5b shows spikes similar to those observed in local simula¬ 
tions. Sano & Inutsuka (2001) showed by local three-dimensional resistive MHD simulations 
that these spikes are created by current dissipation in channel hows which appear in the non¬ 
linear stage of the MRI. Although the BSS magnetic helds which appear in our simulation 
are not exactly the same as channel hows in a disk initially threaded by vertical magnetic 
helds, the BSS helds have current layers where the magnetic held lines change their direc¬ 
tion. Magnetic reconnection can take place in such current layers. We call the region as BSS 
channel. 

Figure 6 shows the distribution of density, pressure, radial velocity and specihc angular 
momentum near the equatorial plane averaged azimuthally, vertically (0 < < 0.3) and 
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in time. Solid curves and dashed curves show average in 29000 < t < 31000 and 28000 < 
t < 29000, respectively. Around zu ~ 10, The radial structure changes from a dense torus 
to an accreting disk (hgure 6a). The equatorial density increases with time because the 
torus flattens. Figure 6b shows the averaged gas pressure and the magnetic pressure. In the 
innermost region (2 < w < 3), magnetic pressure is about 30% of the gas pressure. Figure 
6 c shows the distribution of the radial velocity. In the inflowing region {zu < 10), the radial 
velocity is roughly proportional to zu~^ and exceeds 0.1c around the radius of marginally 
stable orbit {zu ~ 3). The radial velocity in the outer region [zu > 10) is much smaller 
than the sound speed. Figure 6d shows the specihc angular momentum averaged in the 
vertical direction, azimuthal direction and in time. The dotted curve shows the Keplerian 
angular momentum distribution. Although the distribution of initial angular momentum is 
uniform, the angular momentum distribution becomes nearly Keplerian in the radial range 
3 < ci7 < 100. 

Figure 7 shows the time development of the density averaged azimuthally and vertically 
(0 < 2 : < 0.3). Gray scale shows the logarithmic density whose range is —0.5 (blue) < log p < 
0 (pink). The interface between the circularly rotating dense disk and radially infalling flow 
(blue region) moves inward from cjj ~ 30 at f ~ 24000 to ci7 ~ 3 at f ~ 26000. The 
radius of the interface oscillates around ci7 ~ 10. The infall of dense matter from the inner 
region of the torus {zo ~ 20) takes place intermittently at f ~ 26000 and t ~ 29000. The 
interval of the infall is nearly equal to the rotation period of the dense torus (to ~ 3079). 
The intermittency arises from the oscillation of the torus excited by the angular momentum 
exchange between the infalling matter and the torus material; the dense torus is kicked by 
the infalling blob due to the angular momentum conservation. Note that in hgure 7 the 
location of the density maximum of the torus is pushed outward following the infall of dense 
blobs. Short time scale oscillations in the inner region are overlaid to this global oscillation. 
Figure 8 shows the time development of the azimuthally averaged radial velocity relative to 
the mean radial velocity approximated as ~ max (3ci7“^, 0.003) at equatorial plane. 

Bright region and dark region show the region where > 0 and < 0, 

respectively. The innermost region (2 < ccj < 10) shows quasi-periodic oscillations with 
frequency / ~ fi;max/27r ~ 114O(M/M0)“^[Hz] where Kmax is the maximum of epicyclic 
frequency as already pointed out by one-dimensional simulations (e.g., Matsumoto, Kato & 
Houma 1988; Kato, Fukue & Mineshige 1998). Signihcant inhow occurs around t = 25000, 
27000, 28000, 29500 and 31000. 

Figure 9 shows the position of the Lagrange test particles uniformly distributed in 
the equatorial plane inside ru < 18 at f = 30319. Spiral arms which appear in hgure 9 
correspond to the dense arms in the density distribution (see hgure 3). The diamond shaped 
symbols trace the same test particles. The motion of the particles show the development 



of non-axisymmetric MRI. The particle jl4, for example, has lower angular momentum than 
other particles and spirally infalls. The infalling matter deforms the magnetic helds into the 
BSS shape. We can also see the train of spirally infalling particles around a; = 3, ?/ = 0 at 
t = 30550 and t = 30630. Figure 8 indicates that accretion to the black hole is not time steady 
but intermittent. Figure 8 shows that evacuated region is created around a; = 5,?/ = —5 at 
t = 30590. The sharp ridge of the density distribution around a; = 0, y = —7 in hgure 3b 
indicates the interaction between the outgoing wave and the infalling matter. In order to 
show the propagation of outgoing wave, we show in figure 10 the distribution of density. At 
t = 30570, dense region (colored in orange) appears and propagates outward. As denoted 
by arrows, density enhancement develops in the evacuated region and propagates in the 
direction opposite to the rotation of the disk. In the next subsection, we would like to show 
that the motions of these density enhancements are due to magnetic reconnection. 


3.3. Magnetic Reconnection in the Plunging Region 

Figure 11 shows the time evolution of the current density distribution and magnetic 
held lines in the innermost region (—10 < a; < 10, —5 < y < 5). In the outer region where 
magnetic held lines show turbulent structure (hgure 3b), isocontours of current density have 
fractal like distribution (Kawaguchi et ah 2000). On the other hand, long current sheets 
are formed in the inner infalling region {w < 8) where magnetic held lines show BSS shape 
(hgure 3d). Magnetic energy is accumulated in these current sheets. Red regions in hgure 
11 are active regions where the current density is high. Active region ”A” is a narrow peak 
in the current distribution (current sheet). White arrows denoted by ”A” trace the motion 
of this current sheet. Active region ”B” has a cusp in the current distribution. White 
arrows denoted by ”B” indicate the tip of the reconnected magnetic held lines behind the 
reconnection point. At t = 30590, the current layer in active region ”B” coincides with the 
density enhancement shown in hgure 10b. The lower panels in hgure 11 show that the cusp¬ 
shaped current layer ”B” propagates in the direction relatively opposite to the direction of 
disk rotation and that current dissipation takes place. These are the magnetic haring sites 
where magnetic energy is released by magnetic reconnection. The double peak in the current 
distribution which splits from the cusp is a characteristic feature of slow shocks associated 
with magnetic reconnection (see Priest 1982). 

Figure 12 shows the time evolution of magnetic held lines. Closed magnetic loops 
(magnetic islands) formed around the cusp of magnetic held lines indicate that magnetic 
reconnection is taking place in the current sheet. Figure 12 indicates that magnetic recon¬ 
nection in an accretion disk has two types; (a) magnetic reconnection taking place inside 
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an elongating magnetic loop and (b) magnetic reconnection driven by the interaction of two 
magnetic loops. In figure 13, we schematically show the magnetic reconnections in solar 
corona and in accretion disks. Arcade type solar flares (e.g., Tsuneta et al. 1992) which 
accompany loop eruption belong to the loop elongation type. Solar microflares driven by 
loop interaction produces X-ray brightning and X-ray jets (see e.g., Shibata et al. 1992). 
The loop elongation type reconnection can be seen in active region ”A” at t = 30610 in figure 
11 and figure 12c. Backward propagating reconnection jet appears as indicated by arrows 
in hgure 10c. The active region ”B” in hgure 11 subjects to both type of reconnection. At 
t = 30590, the active region ”B” collides with the preceding magnetic loop (hgure 12b) and 
produces intense current layers. Magnetic reconnection taking place in this current layer re¬ 
leased more energy than active region ”A” because the current layer ”B” has larger volume. 
Although loop interaction releases less energy than arcade eruptions in solar corona, loop 
interaction can release more energy in accretion disks because rotational energy of the disk 
is converted to the magnetic energy by such interactions. Figure 14 shows the time evolution 
of temperature. Magnetic energy is deposited into the thermal energy and heats the plasma 
around the reconnection site. Filamentary shaped hot region which coincide with current 
sheets trace the slow shocks associated with magnetic reconnection. 

Figure 15a shows the time evolution of the volume integrated Joule heating rate Ej = 
J Tjj'^dV in 2 < CC7 < 6 and 0 < < 10. Figure 15b shows the time evolution of J j‘^dV. 

Figure 15c shows the magnetic energy Em = / /8n)dV. Figure 15d shows accretion rate 
at CtJ = 2.5. Arrows indicate t = 30590 and t = 30630. The largest flare occurs around 
this time interval. Since we turn on anomalous resistivity only when j / p> {j / p)ci magnetic 
energy is released after the accretion rate decreases and the innermost region is rarehed. The 
Joule dissipation rate pp tends to be zero when the matter density of the BSS channel is 
high. This is the reason why Ej has maximum after the peak of J j'^dV. Magnetic energy 
decreases between t = 30590 and t = 30630 as expected from the magnetic reconnection 
model. Only a small fraction (< 1%) of the released magnetic energy goes into the Joule 
heating because p < 10“^ in this simulation. The released magnetic energy is converted to 
the thermal energy by slow shock (see hgure 14). The volume integrated current density 
correlates well with the mass accretion rate. This is because magnetic helds are stretched 
and amplihed by mass accretion. In other words, parts of the gravitational energy of the 
accreting matter are stored in the current sheet and released by magnetic reconnection. 

Next, let us compare the X-ray light curve obtained by numerical simulations with 
observations of X-ray shots. Figure 16 shows the time variation of the X-ray luminosity 
computed by the optically thin bremsstrahlung luminosity E^ = f p^T^B (dashed 

curve) and Joule heating rate Ej = j pj'^dV (solid curve) where the integration is over the 
innermost region {2 < w < < z < 10). We introduce cutoff function /(T) where 
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f(T) = 1 when 0.01 < T < 0.1 and otherwise f(T) = 0. The vertical scale of hgure 16 
is arbitrary because and Ej depend on the unit density pb and the electric resistivity 
Pq (see equation (5)). As the matter in dense spiral arms accrete to the innermost region, 
X-ray luminosity gradually increases and has a peak around t ~ 30100. The time scale 
of the increase of X-ray luminosity is the accretion time scale of dense blobs and typically 
~ lO^r^/c ~ 0.01M/Mq sec in our simulation (see figure 7). As the dense blob is swallowed 
into the black hole, X-ray luminosity decreases. Subsequently, magnetic reconnection takes 
place in the rarefied region and releases the magnetic energy. The soft X-ray excess precedes 
the magnetic reconnection about SOOr^/c ~ 5M/M0ms. 


4. Summary and Discussion 

We have shown that in the innermost plunging region of black hole accretion disks, 
growth of non-axisymmetric MRI creates bisymmetric spiral (BSS) magnetic fields. Mass 
accretion proceeds along these spiral channels. As a result, density distribution tends to 
show one armed spiral structure. 

A current sheet is created inside the BSS channel because the magnetic held changes its 
direction inside the channel. The current sheet is built up by converting the gravitational 
energy of the accreting gas into magnetic energy. In other words, magnetic energy is accumu¬ 
lated by mass accretion. When the BSS channel is rarehed after dense blobs infall, magnetic 
reconnection takes place and converts the magnetic energy into heats. This mechanism is 
analogous to the accumulation of magnetic energy by mass motion preceding solar hares and 
protostellar hares (e.g., Hayashi et al. 1996). Magnetic reconnection can also take place by 
interaction of two long current sheets. These X-ray hares in the plunging region can be the 
origin of large amplitude X-ray time variabilities characteristic of black hole candidates. 

Let us compare our numerical results with X-ray observations of Cyg X-1 (Negoro 
et al. 1995). In X-ray hard states. X-ray hux from Cyg X-1 shows largest shots whose 
interval is typically several seconds (Negoro & Mineshige 2002). In our simulation, the 
largest shot occurred at O.SM/M© sec after the initial torus is set up at w ~ 50. The 
magnetic reconnection produces outgoing waves and cusp-shaped current sheets. We have 
shown that the mass accretion from the circularly rotating dense torus is not time steady but 
intermittent. The interval of mass feeding from the torus is the rotation period of the torus 
(hgure 7a). In our simulation, this time scale corresponds to to ~ SOOOr^/c ~ 0.03M/Mq 
sec. If the torus locates at = bOOr^, this interval ~ 0.9M/Mq sec approaches the 
interval of observed shots if the mass of the black hole is IOMq. Once dense blobs accrete to 
the innermost region, the prohle of the shot is essentially determined by the local processes 
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in the innermost region. Thus the prohle of one shot may not depend on the location of the 
center of the torus. However, the interval of the shot depends on Smaller X-ray shots 
and fluctuations can be created by magnetic reconnections ubiquitous in accretion disks. 
Observations of Cyg X-1 indicate that the averaged shot prohle is time symmetric (Negoro 
et ah 2001) and X-rays become hard at the peak of the shot. Figure 17 schematically 
shows the mechanism of X-ray shots. The upper left panel shows the accretion stage when 
dense blobs fall into the innermost region. As the density in the innermost region increases, 
the soft X-ray luminosity increases. Current sheets are formed inside the BSS channel and 
magnetic energy is accumulated. After the blobs are swallowed into the black hole, the 
current sheet in the BSS channel is rarehed. The soft X-ray luminosity begins to decrease 
because X-ray emitting gas is depleted. Largest magnetic reconnection takes place in the 
rarehed current sheet by loop elongation or by loop interaction. Since the released magnetic 
energy is converted to heat. X-ray luminosity will suddenly increase and decays exponentially. 
In order to compute the light curve after magnetic reconnection, we will have to include the 
effects of heat conduction along the magnetic loop which may evaporate the disk material 
(see the numerical simulation of solar hares by Yokoyama & Shibata 2001). We would like 
to report the results of such simulations in future. 

The BSS magnetic helds in the innermost region of black hole accretion disks contribute 
to the angular momentum transport of the disk material. In this region, angular momen¬ 
tum transport due to this ordered spiral magnetic helds is more efficient than the angular 
momentum transport by turbulent magnetic helds. 

The dynamical ehects of magnetic helds in the plunging region of black hole accretion 
hows have been investigated analytically by Krolik (1999), Gammie (1999), Agol & Krolik 
(1998, 2000) and simulated by Hawley & Krolik (2001, 2002) and Krolik & Hawley (2002). 
Our numerical results support their conclusion that the ratio of stress to pressure (a) has 
systematic gradient with radius and has larger values well inside the plunging region. 

In the innermost region of the disk, the accretion how is dominated by radial advection. 
The radial structure of the innermost region {w < lOr^) of accretion how obtained by 
three-dimensional MHD simulations approaches the global transonic solution of optically 
thin disks with the viscous parameter a ~ 0.1 (e.g., Narayan et al. 1997). The direct 
numerical simulations such as the one we presented in this paper have the potential to do 
much more than seek agreement with a-models. They can predict the time varitation of 
black hole accretion hows without assuming the phenomenological a-parameter. In order to 
conhrm the applicability of numerically obtained accretion hows to black hole candidates, it 
is essential to compute the X-ray spectrum from numerical results and compare them with 
observations. We would like to report the results of such analysis in subsequent papers. 



- 12 - 


The authors thank Drs. S. Mineshige, K. Shibata, H. Negoro, T. Sano, and R. Narayan 
for discussion. The authors thank the Yukawa Institute for Theoretical Physics at Kyoto 
University. Discussions during the YITP workshop YITP-W-01-17 on ’’Black Holes, Grav¬ 
itational Lens, and Gamma-Ray Bursts” were useful to complete this work. Numerical 
computations were carried out by VPP5000 at NAOJ. This work was supported in part by 
Research Fellowships of the Japan Society for the Promotion of Science for Young Scientists 
(13011092, MM), and by AGT-JST (RM) of Japan Science and Technology corporation. 



- 13 - 

REFERENCES 

Abramowicz, M. A., Czerny, B., Lasota, J. P., & Szuszkiewicz, E. 1988, ApJ, 332, 646 
Agol, E., & Krolik, J., H. 1998, ApJ, 507, 304 
Agol, E., & Krolik, J., H. 2000, ApJ, 528, 161 
Armitage, P.J. 1998, ApJ, 501, L189 

Armitage, P. J., Reynolds, C. S., & Chiang, J. 2001, ApJ, 548, 868 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 

Brandenburg, A., Nordlund, A., Stein, R. F., & Torkelsson, U. 1995, ApJ, 446, 741 
Gammie, C. F. 1999, ApJ, 522, L57 

Hayashi, M. R., Shibata, K., & Matsumoto, R. 1996, ApJ, 468, L37 

Hawley, J. F., Gammie, C. F., & Balbus, S. A. 1995, ApJ, 440, 742 

Hawley, J. F. 2000, ApJ, 528, 462 

Hawley, J. F., & Krolik, J. H. 2001, ApJ, 548, 348 

Hawley, J. F., Balbus, S. A., & Stone, J. M. 2001, ApJ, 554, L49 

Hawley, J. F. 2001, ApJ, 554, 534 

Hawley, J. F., & Krolik, J. H. 2002, ApJ, 566, 164 

Hawley, J. F., & Balbus, S. A. 2002, ApJ, 573, 738 

Kato, S., Fukue, J., & Mineshige, S. 1998, Black-Hole Accretion Disks, (Kyoto: Kyoto 
University Press) 

Kawaguchi, T., Mineshige, S., Machida, M., Matsumoto, R., & Shibata, K. 2000, PASJ, 52, 
LI 

Krolik, J., H. 1999, ApJ, 515, L73 

Krolik, J. H., & Hawley, J. F. 2002, ApJ, 573, 754 

Machida, M., Hayashi, M. R., & Matsumoto, R. 2000, ApJ, 532, L67 

Machida, M., Matsumoto, R., & Mineshige, S. 2001, PASJ, 53, LI 



- 14- 

Manmoto, T., Takeuchi, M., Mineshige, S., Matsumoto, R., & Negoro, H. 1996, ApJ, 464, 
L135 

Matsumoto, R., Kato, S., Fukue, J., & Okazaki, A. T. 1984, PASJ, 36, 71 

Matsumoto, R., Kato, S., & Fukue, J. 1985, in Theoretical Aspects on Structure, Activity, 
and Evolution of Galaxies III, eds S. Aoki, M. lye, Y. Yoshii (Tokyo: Tokyo Astr. 
Obs. Pubh), p 102 

Matsumoto, R., Kato, S., & Honma, F. 1988, in Physics of Neutron Stars and Black Holes 
ed. Y. Tanaka (Tokyo: Universal Academy Press), p.l55 

Matsumoto, R., & Tajima, T. 1995, ApJ, 445, 767 

Matsumoto, R. 1999, in Numerical Astrophysics, ed. S. M. Miyama, K. Tomisaka, & T. 
Hanawa (Amsterdam: Kluwer Academic Publishers), p.l95 

Muchotrzeb, B & Paczyhski, B. 1982, Acta Astronomica, 32, 1 

Narayan, R., Kato, S., & Honma, F. 1997, ApJ, 476, 49 

Negoro, H., Kitamoto, S., Takeuchi, M., & Mineshige, S. 1995, ApJ, 452, L49 

Negoro, H., Kitamoto, S., & Mineshige, S. 2001, ApJ, 554, 528 

Negoro, H., & Mineshige, S. 2002, PASJ, 54, L69 

Okada, R., Fukue, J., & Matsumoto, R. 1989, PASJ, 41, 133 

Paczyhski, B., & Witta, P. J. 1980, A&A, 88, 23 

Priest, E., R. 1982, Solar Magneto-Hydrodynamics (Boston: D.Reidel Publising Company), 
p.ll4 

Reynolds, C., S., & Armitage, P. J. 2001, ApJ, 561, L81 
Sano, T., & Inutsuka, S. 2001, ApJ, 561, L179 

Shibata, K., Ishido, Y., Acton, L. W., Strong, K. T., Hirayama, T., Uchida, Y., McAllister, 
A. H., Matsumoto, R., Tsuneta, S., Shimizu, T., Hara, H., Sakurai, T., Ichimoto, K., 
Nishino, Y., & Ogawara, Y. 1992, PASJ, 44, L173 

Tajima, T., & Gilden, D., 1987, ApJ, 320, 741 



- 15 - 

Tsuneta, S., Kara, H., Shimizu, T., Acton, L. W., Strong, K. T., Hudson, H. S., & Ogawara, 
Y. 1992, PASJ, 44, L63 

Yokoyama, T., & Shibata, K. 1994 ApJ, 436, L197 
Yokoyama, T., & Shibata, K. 2001 ApJ, 549, 1160 


This preprint was prepared with the AAS MpK macros v5.0. 



16 


Fig. 1.— (a) Initial model of the torus, (b) The isosurface of density (p = 0.4) at t = 30590. 
Numbers show the (x, p, z) coordinate. 

Fig. 2.— (a) Azimuthally averaged density distribution logp (color scale) and the poloidal 
momentum vectors (arrows) at t = 30590. Long arrows are plotted in dark color, (b) 
Azimuthally averaged density distribution at t = 0, t = 21502 and t = 30590. 

Fig. 3.— The equatorial density distribution (color scale) and magnetic held lines projected 
onto the equatorial plane (grey curves) at t = 30590. (a), (b) Global structure inside GOr^. 
(c), (d) Inner region inside lOr^. Color scale of logarithmic density logp is shown at the 
bottom of each panel. 

Fig. 4.— Three-dimensional magnetic held structure at t = 30590. White curves show 
magnetic held lines, blue surface show the isosurface of the density (p = 0.4), and yellow 
plane shows the contour of the equatorial density. The size of the box is —10 < a; < 10, 
— 10 < p < 10 and —10 < 2 : < 10 . 

Fig. 5.— (a) The time developments of the volume averaged plasma (3 = (P)/ /Stt) . 
(b) The time developments of the ratio of the volume averaged Maxwell stress to pressure, 
as = —{B^B^/An)/(P). (c) The time developments of the accretion rate at zu = 2.5rg. The 
solid curves show the average in the inner region 4 < ccj < 10 and 0 < < 1. The dotted 

curves show the average in the outer region of the disk, 20 < ciJ < 40, and 0 < < 3. 

Fig. 6 .— Radial prohles of density, gas pressure, magnetic pressure, radial velocity and 
specihc angular momentum averaged in azimuthal direction, vertical direction (0 < < 0.3) 

and in time. Solid curves show quantities averaged in 29000 < t < 31000 and dashed curves 
show quantities averaged in 28000 < t < 29000. The dotted curve in (d) shows the Keplerian 
angular momentum distribution. 

Fig. 7.— Time development of the azimuthally averaged density at the equatorial plane. 
The range of the color scale is —0.5 (blue) < logp < 0 (pink). 


Fig. 8.— Time development of the azimuthally averaged radial velocity at the equatorial 
plane Av^ = where (v^) = — max (3ci7“^, 0.003). The range of the gray scale is 

—0.07 (dark) < Av^ < 0.07 (bright). 
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Fig. 9.— Positions of Lagrange test particles on the eqnatorial plane. Large symbols (jjl ~ 
jl6) trace the same particles. 


Fig. 10.— Density distribntion in the region 2 < zu < 10. Color scale in each panel shows 
the scale of logp. Numbers at the upper right corner of each panel show time. Arrows 
indicate the reconnection jet. 


Fig. 11.— Time evolution of the current density (color scale) and magnetic held lines (white 
curves). Red area show active regions where the current density is high. White arrows ”A” 
trace the motion of a current sheet. White arrows ”B” indicate the tip of the reconnected 
magnetic held lines. The size of each box is —10 < a; < 10 and —5 < y < 5. Numbers show 
time. 


Fig. 12.— Magnetic held lines projected onto the equatorial plane. Numbers show time. 
The size of the box is —5 < a; < 5 and —5 < p < 5. 


Fig. 13.— A schematic picture of magnetic reconnection in solar hares and in accretion 
disks. 

Fig. 14.— Distribution of logarithmic equatorial temperature. The size of each box is 
2 < w < 10. 

Fig. 15.— Time evolution of (a) the volume integrated Joule heating rate, (b) the volume 
integrated current density, (c) the magnetic energy, and (d) the accretion rate at the zu = 2.5. 
Arrows indicate t = 30590 and t = 30630. 


Fig. 16.— Time variation of the X-ray luminosity computed by the simulation. Solid 
curve shows the Joule heating rate. Dashed curve shows the soft X-ray luminosity by 
bremsstrahlung. 


Fig. 17.— A schematic picture of the mechanism of X-ray shots. Soft X-ray excess precedes 
the hard X-ray hare. 
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